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Abstract 



Wc give a complete characterization of the two-state anti-ferromagnetic spin systems which 
are of strong spatial mixing on general graphs. We show that a two-state anti-ferromagnetic 
spin system is of strong spatial mixing on all graphs of maximum degree at most A if and only 
if the system has a unique Gibbs measure on infinite regular trees of degree up to A, where 
A can be either bounded or unbounded. As a consequence, there exists an FPTAS for the 
partition function of a two-state anti-ferromagnetic spin system on graphs of maximum degree 
at most A when the uniqueness condition is satisfied on infinite regular trees of degree up to 
A. In particular, an FPTAS exists for arbitrary graphs if the uniqueness is satisfied on all 
infinite regular trees. This covers as special cases all previous algorithmic results for two-state 
anti-ferromagnetic systems on general-structure graphs. 

Combining with the FPRAS for two-state ferromagnetic spin systems of Jcrrum-Sinclair and 
Goldberg- Jerrum-Paterson, and the very recent hardness results of Sly-Sun and independently of 
Galanis-Stefankovic-Vigoda, this gives a complete classification, except at the phase transition 
boundary, of the approximability of all two-state spin systems, on either degree-bounded families 
of graphs or family of all graphs. 
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1 Introduction 



Spin systems are well studied in the areas of Statistical Physics, Applied Probability and Computer 
Science as a general framework to capture the essence of how local interactions and constrains affect 
the macroscopic properties of particle systems. A system is usually described by a graph, with each 
vertex in one of a fixed number of states called spins, and edges specifying the neighborhood relation 
of the system. 

Let G(y, E) be a graph and q be the number of spin states. A configuration of the system is one 
of the gl^l possible assignments a : V ^ [q\. Each configuration a has an energy H{a) as a sum over 
all edges and vertices, such that the contribution of an edge {u, v) € E is determined by a symmetric 
function of the spin states cr(u) and cr{v), and the contribution of a vertex v £ V is determined 
by a function of its spin state cr{v). The weight of a configuration a is w{a) = exp(— ■^^^), where 
T is the temperature. We focus on the two-state spin systems. Up to normalization, a two-state 
spin system is fully captured by three parameters (/3,7, A), where /3 and 7 determine the symmetric 
function for edge contribution and A, also known as the external field, determines the function for 
vertex contribution. The Gibbs measure is a natural probability distribution over all configurations 
such that the probability of a configuration a is where the normalizing factor Z = ^w(cj) 

is called the partition function. The partition function encodes rich information regarding the 
macroscopic behavior of the spin system. However, for almost all nontrivial settings it is ^^P-hard 
to compute the precise value of partition functions. 

One of the most important properties of spin systems is the correlation decay, which says that 
the correlation between the marginal Gibbs distributions of two vertices decays rapidly with respect 
to the distance between the two vertices. This property is also called weak spatial mixing [25]. Of 
greater algorithmic significance is the strong spatial mixing, which says that the correlations decay 
in the presence of arbitrary fixed spins at other vertices. For two-state spin systems, the strong 
spatial mixing may imply efficient approximation algorithms for the partition function. It is then 
an important question to characterize the systems which exhibit strong spatial mixing on arbitrary 
instances of graphs by the parameters of the systems. 

A two-state spin system is ferromagnetic if adjacent vertices favor agreement of spins, and is 
anti-ferromagnetic if otherwise. For all two-state ferromagnetic spin systems, the partition functions 
can be efficiently approximated jl0pi2j . In the anti- ferromagnetic region, the correlation decay plays 
a central role in the approximability of partition functions. It is believed (see [18| ) that for such 
models the approximability of the partition function is characterized by the uniqueness of Gibbs 
measure on infinite regular trees, which is equivalent to the weak spatial mixing on the infinite 
regular trees. This condition is called the uniqueness condition. 

In a seminal work |25] , Weitz shows that for the hardcore model strong spatial mixing is charac- 
terized by the uniqueness condition, and in a recent work of Sinclair, Srivastava, and Thurley |20] 
the same characterization is proved for the anti-ferromagnetic Ising model. Both models are im- 
portant special two-state anti-ferromagnetic spin systems. On the hardness side, it is proved for 
the hardcore model by Sly [22] and Galanis et al. [6], and very recently for the general two-state 
anti-ferromagnetic spin systems by Sly and Sun [23] and independently by Galanis, Stefankovic, 
and Vigoda [7| that violating the uniqueness condition implies the inapproximability of partition 
functions. Two questions remain open for our complete understanding of the correlation decay 
and computation in two-state spin systems: the characterizations of the strong spatial mixing and 
approximability of general two-state spin systems. 
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1.1 Our results 



We characterize the two-state anti-ferromagnetic spin systems which exhibit strong spatial mixing 
on general degree-bounded graphs or arbitrary graphs by the uniqueness of Gibbs measure on 
infinite regular trees. 

Theorem 1. For any finite A>2 or A = oo, a two-state anti-ferromagnetic spin system is of 
strong spatial mixing on graphs of maximum degree at most A if and only if the system exhibits 
uniqueness on infinite d-regular trees for all d < A. 

Due to Weitz's self-avoiding walk tree construction [25], the strong spatial mixing of a two-state 
spin system on degree-bounded graphs immediately implies an FPTAS for the partition function. 
Indeed, we show an even stronger notion of correlation decay introduced in a previous work |15j . 
namely the computationally efficient correlation decay, which gives FPTAS not only for the degree- 
bounded graphs but also for arbitrary graphs with unbounded degrees, when the corresponding 
uniqueness condition is satisfied. 

Theorem 2. For any finite A > 3 or A = oo, there exists an FPTAS for the partition function 
of the two-state anti- ferromagnetic spin system on graphs of maximum degree at most A if for all 
d < A the system parameters lie in the interior of the uniqueness region of infinite d-regular tree. 

In the above two theorems, the A = oo case represents the graphs of unbounded degree. 

Due to a very recent hardness results of Sly and Sun [23] for general two-state spin systems and 
an independent result of Galanis, Stefankovic and Vigoda [7] for a less general setting, violating 
the uniqueness condition implies inapproximability of the partition function. 

Theorem 3 (Due to [23] and [7]). For any finite A > 3 or A = oo, unless NP = RP there does 
not exist an FPRAS for the partition function of the two-state anti-ferromagnetic spin system on 
graphs of maximum degree at most A if for some d < A the system parameters lie in the interior 
of the non-uniqueness region of infinite d-regular tree. 

The original theorem in [23] holds for d-regular graphs with fixed d, which immediately implies 
the hardness for degree-bounded graphs or arbitrary graphs. And the hardness condition in both 
|23j and [7] requires the non-uniqueness on a d-regular tree with d > 3. But the uniqueness on the 
infinite 2-regular tree (i.e. the infinite path) always holds for any two-state anti-ferromagnetic spin 
system, thus the condition in Theorem [3] suffices. 

For graphs of maximum degree 2 or less, the partition function can be computed exactly in 
polynomial time. Thus Theorem [2] and Theorem [3] together with the FPRAS for two-state ferro- 
magnetic spin system [THllIl] give an almost complete (except at the phase transition boundary) 
classification of the approximability of the partition function of all two-state spin systems, on either 
degree-bounded families of graphs or family of all graphs. 

1.1.1 Regularity and monotonicity 

In Statistical Physics, the correlation decay is usually studied on regular graphs or even structurally 
symmetric graphs (e.g. Bethe lattice). Although for hardness results considering only regular graphs 
will strengthen the lower bounds, from algorithmic perspective it is more general to consider spin 
systems on general graphs. The approximation algorithms in [25] for the hardcore model and [20] 
for the anti- ferromagnetic Ising model are both for general graphs with bounded maximum degrees. 
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As discussed in |20y23] , up to translation of parameters the hardcore and Ising models together are 
complete for general two-state spin systems on d-regular graphs with fixed d. However, this does 
not cover the most general case, namely the general two-state spin system on general graphs, of 
either bounded or unbounded degree. A fundamental reason for this gap is the non-monotonicity 
of general spin systems. 

The well studied hardcore and anti-ferromagnetic Ising models (along with all two-state spin 
systems with /3,7 < 1) are both monotone spin systems, in the sense that the uniqueness on 
infinite d-regular tree implies the uniqueness on all infinite regular trees with smaller degrees. This 
monotonicity does not necessarily hold in general two-state spin systems. 

In |25j . Weitz established the following implication in the hardcore model: 

Claim (Theorem 2.3 in |25]). Strong spatial mixing on a d-regular tree implies the strong spatial 
mixing on all graphs of maximum degree at most d. 

In |25j . Weitz also remarked without proof that this implication holds for all two-state spin 
systems (indeed this is rigorously proved for anti- ferromagnetic Ising model in [20]). With this to 
be true, devising approximation algorithms for two-state spin systems on degree-bounded graphs is 
reduced to verifying the strong spatial mixing on d-regular trees. This has been accepted as a fact 
about the two-state spin systems and has become a building block for approximation algorithms 
for such systems (see Theorem 2.8 in [20] and the discussions in [211122] ). It was also raised as a 
conjecture in [21] whether the claim holds for general multi-state spin systems. 

As a byproduct of our analysis (see Section [5|), we find that this well-believed implication 
between the strong spatial mixing on d-regular tree and on graphs of maximum degree at most d 
holds only for monotone spin systems (including the hardcore and anti-ferromagnetic Ising models). 
For general two-state spin systems the worst case for uniqueness as well as strong spatial mixing 
among all degree-bounded graphs is indeed a regular tree, but may no longer be the one of the 
highest degree. A bright side of this complication is that higher degrees may yield much faster 
correlation decay, making possible the FPTAS for graphs with unbounded degrees. 

These new phenomena suggest that the general two-state spin systems have much richer struc- 
ture than the well-studied monotone spin systems such as the hardcore and anti-ferromagnetic 
Ising models. The former approach via the strong spatial mixing on d-regular tree which succeeds 
in monotone spin systems on general graphs and general spin systems on regular graphs, meets 
a barrier when dealing with general spin systems on general graphs. We give a unified approach 
to the correlation decay in general two-state spin systems, through the strong spatial mixing on 
arbitrary trees instead of d-regular trees. This approach was initiated in our previous work [15] 
dealing with graphs of unbounded degrees. In this paper, we devise a unified potential-based anal- 
ysis which adapts to both the irregularity of the arbitrary tree and the non-monotonicity of general 
two-state spin system and give tight correlation decay results for all two-state anti-ferromagnetic 
spin systems on degree-bounded families of graphs and family of all graphs. 

1.1.2 Implications of the main results 

By solving the uniqueness condition we can restate our main results in various threshold forms. 
Theorem [2] covers as special cases all previous algorithmic results for two-state anti- ferromagnetic 
spin systems on general-structure graphs as well as clears up previously uncovered cases. 
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In terms of interactions. We can fix the external field A and discuss the tractable region of 
(/3,7). Since the roles of /3 and 7 are symmetric, we can further fix one of them and discuss the 
tractable range of the other. This formulation was used in [lOy iSj. 

Our main result can be restated as follows: for any A, there is a critical threshold 7c(/3, A, A) 
for the uniqueness on infinite regular trees up to degree A such that if 7c(/3, A, A) < 7 < there is 
an FPTAS for graphs of maximum degree at most A; and in particular, 7c = 7c(/5, A, 00) > 1 is an 
absolute constant such that if 7 G (7c, ^) there is an FPTAS for arbitrary graphs. 

This covers as special cases all algorithmic results in [10] regarding the anti-ferromagnetic 
systems and all results in [15], extends the tractable regions in these previous works, and considers 
the degree-bounded graphs as well. 

In terms of external field. Motivated by the studies of hardcore and anti-ferromagnetic Ising 
models, we can fix (/3, 7) and discuss the tractable range of external field. 

Due to the symmetric role of /3 and 7, we may assume that /3 < 7. Our main results can be 
restated in specific settings as follows: 

• Hard constraints (when /3 = 0): For any A, Ac(7, A) = mini<rf<A is a critical thresh- 
old for the uniqueness on infinite regular trees up to degree A such that if A < Ac (7, A) there 
exists an FPTAS for graphs of maximum degree at most A. 

For 7 < 1, the critical threshold equals Ac(7, A) = ^ — ' There is no external field 

A > satisfying uniqueness on infinite regular trees of unbounded degrees. This is consistent 
with the hardness result for the hardcore model without degree bound [31l22j . One particularly 
interesting special case is when 7 = 1, in which case the model is exactly the hardcore model 
with fugacity parameter A, and Ac(l, A) = ^^^^^^^^ is the critical fugacity of the uniqueness 
threshold. This covers the result of |25| . 

For 7 > 1, in addition to the results for degree-bounded graphs, there exists an absolute 
positive constant A(7) = mind>i which lower bounds A(7, A) for all A, such that if 

A < A (7) there exists an FPTAS for arbitrary graphs. This is quite different from the case 
7 < 1 and is not covered by any previous results. 

• Soft constraints (when f3 > 0): For any A, Ac = Ac(/3,7, A) and Ac = Ac(/3, 7, A) are two 
critical thresholds for the uniqueness on infinite regular trees up to degree A such that if 
A € (0, Ac) U (Ac, 00) there exists an FPTAS for graphs of maximum degree at most A. In 
particular, if -v/^ > the two ranges overlap, thus for any external field A > there exists 
an FPTAS for graphs of maximum degree at most A. 

When /3 = 7, the spin system becomes the anti-ferromagnetic Ising model. And Ac = 1/Ac, 
thus if A G (0,1/Ac) U (Ac, 00), i.e. |logA| > log Ac, there exists an FPTAS for graphs of 
maximum degree at most A. This covers the result of anti- ferromagnetic Ising model in |20j . 

For unbounded maximum degree A, if 7 < 1, there is no external field A > satisfying 
uniqueness on infinite regular trees of unbounded degrees. This is consistent with the hardness 
results for two-state spin systems in [10]. If 7 > 1, Ac = Ac(/3,7) is a critical threshold for 
uniqueness on all infinite regular trees such that if A < Ac(/3,7) there exists an FPTAS for 
arbitrary graphs. 
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1.2 Related works 



The approximation of partition functions of spin systems has been extensively studied [3H5l[9l[TTI- 
[T3 t ll6 | [M] . In a seminal work |12j , Jerrum and Sinclair devised a fully polynomial-time randomized 
approximation scheme (FPRAS) for the ferromagnetic Ising model. Later in [10], the FPRAS 
was extended to all two-state ferromagnetic spin systems by translating the parameters to the 
ferromagnetic Ising model. Also in [10], Goldberg, Jerrum, and Paterson gave an FPRAS and 
inapproximability results for two-state anti-ferromagnetic spin systems on arbitrary graphs. A 
gadget based on random regular bipartite graphs was proposed by Dyer, Frieze, and Jerrum in [3] 
and was also analyzed by Mossel, Weitz, and Wormald in [18] to study the inapproximability 
on degree-bounded graphs. It is widely believed that the transition of approximability of anti- 
ferromagnetic spin systems is captured by the phase transition of uniqueness on infinite trees. This 
was raised openly as a conjecture in [18]. The conjecture was proved by Sly in [22 for the hardcore 
model. This result was improved by Galanis et al. in [6] to a wide range of parameters. Very 
recently, Sly and Sun |23j proved the hardness of all two-state anti-ferromagnetic spin systems of 
non-uniqueness on infinite regular trees. A same result for anti-ferromagnetic Ising model without 
external field was independently proved by Galanis, Stefankovic, and Vigoda in [7]. 

The correlation decay technique developed independently by Weitz [25] and Bandyopadhyay 
and Gamarnik [T] is a powerful tool for devising deterministic fully polynomial-time approximation 
schemes (FPTAS) for partition functions (other important examples include [2l[8]). In [25], Weitz 
introduced the concept of strong spatial mixing and used it to devise FPTAS for the hardcore 
model up to the uniqueness threshold. The other most important two-state anti-ferromagnetic 
spin system, the anti-ferromagnetic Ising model, was studied recently by Sinclair, Srivastava, and 
Thurley in |20j . where a more powerful message-decay method was introduced to analyze the strong 
spatial mixing and give FPTAS up to uniqueness threshold. A powerful technique was developed 
by Restrepo et al. in [19] which makes use of the specific structure of graphs for strong spatial 
mixing. A broader tractable region than the region of uniqueness is achieved on grid lattice by 
exploiting the structure of the graph. In a previous work [15], we gave an FPTAS for two-state 
anti-ferromagnetic spin systems without external field on arbitrary graphs with unbounded degrees, 
up to a continuous relaxation of the uniqueness threshold. The approach used in the current paper 
was initiated in |15j . however the analysis in [15] cannot separate the uniqueness up to certain 
degree, thus fails in dealing with degree-bounded families of graphs. 



2 Definitions and preliminaries 



A two-state spin system is described by a graph G = {V, E). A configuration of the system is one of 
the 2l^l possible assignments a -.V ^ {0, 1} of states to vertices. We also use two colors blue and 
green to denote these two states. The weight of a configuration can be described as a product of 

^0,0 ^0,1 



contributions of individual edges and vertices. Let A 
of a configuration o":!/— 7'{0,l}is given by 



.^1,0 ^1,1 



and h = (6o , &i ) • The weight 



w{a)= n A 



CF{u),a{y) 



{u,v)eE 



11" 



•a{v) ■ 
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The Gibbs measure is a probability distribution over all configurations defined by p{a) = The 
normalization factor Z{G) = ^f^w{a) is called the partition function. 

We can normalize the contributions of a {blue, green} edge and of a green vertex to be 1. So 

that A = ^ ^ for some /?,7 > 0, and b = (&0)^i) = (-^) 1) for some A > 0. Since the roles of 

blue and green are symmetric, we can assume that /5 < 7 without loss of generality. The three 
parameters (/?, 7, A) with < /? < 7 and A > completely specify a two-state spin system. A 
two-state spin system with /? = 7 is an Ising model and a two-state spin system with /3 = 0, 7 = 1 
or symmetrically /3 = l,7 = 0isa hardcore model. 

A two-state spin system is called anti-ferromagnetic if adjacent vertices favor disagreeing spins, 
i.e. if (3^ < 1. Without loss of generality, we focus on the cases that /3 < 7. 

Definition 4. (/3,7, A) is anti-ferromagnetic «/ < /5 < 7, 7 > 0, /37 < 1, and A > 0. 

By the symmetry of /3 and 7 and the triviality of the case /3 = 7 = 0, this definition is complete 
for all nontrivial two-state anti-ferromagnetic systems 

2.1 Correlation decay 

The Gibbs measure defines a marginal distribution of state for each vertex. Let denote the 
probability of vertex v to be colored blue. Let cja S {0, 1}^ be a configuration of A C V. We call 
vertices v £ A fixed vertices, and v ^ A free vertices. We use p^'^ to denote the marginal probability 
of V to be colored blue conditioning on the configuration of A being fixed as cja- 

Definition 5. A spin system on a family of graphs is said to be of strong spatial mixing if for any 
graph G = {V, E) in the family, any v £V,A dV and a^, t\ G {0, 1}^, 

<exp(-17(distK5))), 

where S C A is the subset on which a\ and ta differ, and dist(u, S) is the shortest distance from v 
to any vertex in S. 

The weak spatial mixing can be defined by measuring the decay with respect to dist(u, A) instead 
of dist(u,S'). The spatial mixing property is also called correlation decay in Statistical Physics. 

Let T be a tree rooted by v. Define R^'^ = Py^ /(l—pZ^) to be the ratio between the probabilities 
that the root v is blue and green, when imposing the condition cja (with the convention that 
= 00 when = 1). Suppose that v has d children and Tj is the subtree rooted by the i-th 
child. Due to the independence of subtrees, we have an easy recursion for calculating R^^: 

Let G{V, E) be a graph. Similarly define that = p^^ /{I —p%^)- In contrast to the case of tree, 
there is no easy recursion for calculating R'o'y for a general graph G because of the dependencies 
caused by cycles. In [2S], a construction called the self-avoiding walk (SAW) tree was introduced 
which reduces the computing of marginal distribution in a general graph to that in a tree. Specif- 
ically, given a graph G{y,E) and a vertex v £ V. The SAW tree Tsaw(G')'^') is a tree rooted 
at V with a new vertex set Vsaw (which effectively enumerates all paths originating from v m. G 
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and may include fixed leaves). Moreover, any vertex sets S C A C V are mapped to respective 
'S'sAW C AsAW C VsAW and any configuration € {0, 1}^ is mapped to a ctxsaw ^ 1}'^'='^'^. We 
abuse the notation and write S = Ssaw and ua = c^Asaw ambiguity is caused. Given a graph 
G{V, E), V £ V and S C V, let distG(f, S) be the shortest distance in G from v to any vertex in S. 

Theorem 6 (Theorem 3.1 of Weitz [25]). Let G{V,E) be a graph, v £V, aA£ {0, 1}^ a configu- 
ration on A C V , and S CV. Let T = TsAwiG,v). It holds that the maximum degree ofT equals 
the maximum degree ofG, distG'('y, S) = distr(f,S'), and Rq^^ = . Moreover, any neighborhood 
of V in T can be constructed in time proportional to the size of the neighborhood. 

2.2 The uniqueness condition 

We consider the uniqueness of Gibbs measure on the infinite {d + l)-regular trees T'^, in which the 
recursion is given by fd{x) — A ^^p^^ due to the symmetric structure of T'^. 

Let Xd be the positive fixed point of fd{x), that is, Xd = f{xd)- It is known [TUITT] that the 
two-state anti-ferromagnetic spin system on T*^ undergoes a phase transition at \fd{xd)\ = 1 with 
uniqueness when \f'^{xd)\ < 1. This motivates the following definition. 

Definition 7. Let (/3, 7, A) be anti- ferromagnetic. Let Xd be the positive fixed point of function 



d 



fd{x) = A j • We say that (/3,7, A) is up-to-A unique, if for all integer 1 < d < A, 

^1 ^ \d{l-P^){Pxd + lY-^ ^ d{l - P^)Xd 

{xd + lY+^ (/3x, + l)(xrf + 7) 

In particular, (/3,7, A) is universally unique if it is up-to-00 unique. 

Being up-to-A unique is equivalent to that the system is of weak spatial mixing on infinite 
regular trees up to degree A. The uniqueness condition can be described in various threshold 
forms, which are given in Appendix lAl 

The uniqueness is defined by requiring that |/^(xd)| < 1. The following lemma states that 
\fd{xd)\ is bounded by an absolute constant as long as the uniqueness condition holds. 

Lemma 8. Let {j3, 7, A) be anti-ferromagnetic. If (/5, 7, A) is up-to-A unique then there exists an 
absolute constant c < 1 which depends only on (3, 7, A and A, such that \fd{xd)\ < c for all 
l<d< A. 

Proof. The lemma holds trivially for finite A. It then remains to show that in case of universal 
uniqueness, |/^(xd)| cannot be arbitrarily close to 1 as grows to infinity. If (/?,7, A) is universally 
unique, due to Lemma 121121 we must have 7 > 1. For anti-ferromagnetic 7, A), /? < ^, thus the 

fixed point Xd = X (^)^ < ^, therefore \f'^{xd)\ = ^pf^;^)jl%^) < The lemma follows. □ 

3 The strong spatial mixing on general graphs 

In this section we prove Theorem [H The necessity of the uniqueness condition is trivial since strong 
spatial mixing on general graphs implies weak spatial mixing on regular trees. It then remains to 
prove the following theorem. 
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Theorem 9. Let (/3,7, A) be anti-ferromagnetic. For any finite A > 2 or A = oo, i/(/3,7,A) is 
up-to-A unique, then the two-state spin system of parameters 7, A) is of strong spatial mixing 
on graphs of maximum degree at most A. 

Our approach is to prove the strong spatial mixing on arbitrary trees of maximum degree at 
most A, which by the self-avoiding walk tree construction implies the theorem. Note that unlike 
in [25] and [20], we analyze the decay on arbitrary tree instead of regular tree. This is because for 
general two-state anti-ferromagnetic spin systems the worst case for strong spatial mixing among 
all graphs of maximum degree at most d may no longer be the d-regular tree. We will explain this 
in detail in Section [5l 

Given any graph G{V,E) of maximum degree at most A, any configuration G {0,1}^ on 
A C V and any S C A, fixing an arbitrary vertex v £ V, hy Theorem [6l a self-avoiding walk 
tree T = TsAwiG,v) can be constructed such that the maximum degree of T is bounded by 



A, distGlf,^) = distT{v,S) and i?^*^ = R^''. Recah that i?^^^, = ^^^/(l -p^'^) thus = 
^'G,v/i^ + ^(tv)- '^A and TA, 



This motivates the following definition. 

Definition 10. Let T be a tree rooted by vertex v, t/^ G {0, 1}^ be a configuration on A C V and 
S" C A 6e a vertex set. Define Ry and 6y as that Ry < RT^ ^ Rv + for all a\ G {0, 1}^ which 
differ from t\ only on S. 

It is then sufficient to prove Theorem [9] by constructing such R^, and 6^, for T = Tsaw(G', v) and 
showing that 5„ < exp(— ri(dist(i), S*))). 

Let T be a tree rooted by v, who has d children vi, . . . ,Vd, and Tj be the subtree rooted by Vi. 
It holds that 



The lower and upper bounds R^, and R^ + 5„ can be recursively constructed as follows. The base 
cases are: (1) v € S, in which case R^ = and 5^ = 00; and (2) v £ A \ S, i.e. v is fixed to be 
the same color in all a\, in which case 6^=0 and R^, = 00 (or = 0) ii v is fixed to be blue (or 
green). For v ^ A, since f{Ri, ■ ■ ■ , Rd) is monotonically decreasing for anti-ferromagnetic (/3,7, A), 



where Ry. and Ry. -\- 6y- are the corresponding lower and upper bounds for R^^ , 1 < i < d. In 
particular, when d = the empty product equals 1 by convention, thus Ry = Ry -\- 5y = A, which 
is consistent with the case that is a free vertex having no children. 

By the monotonicity of f{Ri, . . . ,Rd), it is easy to check that the Ry and Ry + 5y constructed 
above satisfy the requirement of Definition [101 Our goal is to show that 6y decays exponentially 
in depth of recursion when the uniqueness holds. A straightforward approach is trying to prove 




Rrt-} L^rp . 




Ry / {Rvi ~l~ 5 • • • ) -^V^ ~l~ ^V^ ) 1 

Rv ~\~ — fi^Vl ) • • • ) ) ) 
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that 5 contracts at a constant rate in each recursion step. But this does not have to be true to 
guarantee the exponential decay. Indeed there are cases that the error does not decay in single 
steps but decay in a long run. To overcome this, we use a potential $ to amortize the contraction 
and show that 5 ■ $ contracts at a constant rate. We choose the potential function to be 



y/R{l3R+l){R + l) 

We are analyzing the decay on an arbitrary tree with irregular degrees. In order to adapt this 
irregularity, the potential function cannot have d as an input, but only caries the information 
regarding the distribution at the current vertex, yet it has to be able to provide correct compensation 
to the step-wise decay at any state of R and for all spin systems satisfying sufficient uniqueness. A 
heuristic procedure which leads us to this good potential function is discussed in Appendix O 
Let '^{R) be a monotone function satisfying that ^'{R) = ^(R)- We define that 

€y = ip{Rv + 5y) - (p{Rv), 

and accordingly, e^^ = ^{Rvi + ^Vi) - '^{Rvi), I < i < d. 
We define a function a{d;xi, ...,X(i) as follows: 



a{d;xi, ...,Xd] 



A 



The following lemma is obtained from applying the Mean Value Theorem. Similar routines were 
previously used in |15lll9j . The proof of the lemma is in Appendix [Bl 

Lemma 11. The foUowings hold for e„. 

1. (relation to 6^;) = 6^ ■ ^{R) for some R £ [Ry,Ry + 6^]. 

2. (absolute bound) Assuming that 7 > 1 or the maximum degree ofT is bounded by a constant, 
if V ^ A then Ry + 6^ = 0(1) and if Vi 5 for all 1 <i < d then e„ = 0(1). 

3. (stepwise contraction) There exist Ri G [i?^,. + 1 < i < such that 

(-V < a{d;Ri,. . . ,Rd) ■ max{ey-}. 

l<i<d 

To prove the strong spatial mixing, we first relate 5y to by Item [T] of Lemma [TTl and then 
apply induction on the depth in T, with Item [2] of Lemma [TT] as basis and Item [3] of Lemma [TT] as 
induction step. We then need to bound the contraction rate a{d; xi, Xd). Note that < 

(1 + ^ < 1 for 2; G [0,00), thus it holds unconditionally for all Xi £ [0,oo), 1 < i < d, that 

a{d;xi, ...,Xd) < d, (2) 



a{d;xi,...,Xd) < d- ^ (3) 
With the uniqueness, the following much tighter contraction bound can be proved. 
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Lemma 12. Let (/3,7,A) be anti- ferromagnetic. If{(3,j,X) is up-to-A unique, then there exists a 
constant a < 1 such that for any integer 1 < d < A and any xi,...,Xd € [0, +oo), \ < i < d, it 
holds that a{d; xi, . . . , Xd) < a. 

This lemma is the technical core of our analysis. It crucially relies on the choice of potential 
function. Before delving into the formal proof of Lemma [T2l we note that Theorem[9]can be implied 
by this lemma. 

Proof of Theorem [9l 

Let T = TsAw{G,v) for a G whose maximum degree is at most A. Then the maximum degree of 
T is at most A, thus the root v has at most A children in T, and every other vertex in T has less 
than A children. We recursively construct Ru, Su and eu for every subtree in T. 

Let t = dist{v,S). By repeatedly applying Item [3] of Lemma \TT\ without loss of generality, 
we have a path uiU2 ■ ■ ■ ut^2 in T with ui = v such that e^^ < a{dj; xi, . . . , x^j) ■ ^Uj+i for j = 
1,2, ... ,t — 3, where dj is the number of children of Uj and Xi G [0, oo), 1 < i < dj. 

Note that di < A, and dj < A for all other j. Assume that (/3,7, A) is up-to-A unique. If 
A is bounded, then by Lemma [12] there exists a constant a < 1, such that £„. < a ■ e?i ^i for 
2 < j < t — 3, and < di • e„2 < A • 62 due to ([2]); and if A = 00, then by Lemma [T2l e„^. < a • eu^+i 
for all 1 < j < t — 3. In both cases we have = 0{a* ■ eut-2)- 

Due to Item [U of Lemma [TTl 5^ = -7^ = O (—7-^ ■ a*eu. , ) for some R G \Rv,Rv + dv]. We 

then bound ^(R) and eut-i- Due to Item[2]of Lemma [21] the fact that (/3,7, A) is up-to-A unique 
implies that either A is bounded or 7 > 1. Note that v must be free or the theorem is trivial to 
prove, and none of ut_2's children is in 5 because dist(?;,S') = t. Thus by Item [2] of Lemma W\\ 



As discussed in the beginning of this section, this proves Theorem [9l 
Proof of Lemma 1121 

The rest of this section is devoted to the proof of Lemma 1121 Given that (/3,7, A) is up-to-A 
unique, there exists an absolute constant a < 1 such that a(d; x\,. . . , Xd) < a for any 1 < d < A 
and xi, . . . , Xd > 0. 

We define the symmetric version of a{d; xi, . . . , Xd). 



ad{x) = 




The following lemma shows that the symmetric case dominates the maximum of a{d;xi, ...,Xd) by 
using the inequalities of Cauchy-Schwarz and arithmetic and geometric means. 

Lemma 13. Let {/3,-y,X) be anti- ferromagnetic. Then for any integer d and any xi,...,Xd S 

[0,-|-oo), there exists an x £ [0, -|-oo) such that a{d; xi, Xd) < a{d,x). 



e, 




n{i). 

t6^, = 0{a^). 
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Proof. Let Zi = ^^p^- Then Zi G (/3, ^] and Xi = Express a{d;xi, ..■,Xd) in terms of zf. 



a{d;xi,...,Xd) = r r • 2^(^i -'y)Hzi-P)^ 



1 



d _\2 



Due to Cauchy-Schwarz inequality, 



d / 1 \ ^ / 1 ^ 

i=l \ 1=1 / \ i = l / 



Due to the inequahty of arithmetic and geometric means, 



Let z = ( Yli-i Zi j . Then combining the above calculations. 



{Xz^)-2 ■ d{l + - ^z - Pz-')'2 _ I Xzdjz-^ - j){z - ^) 
' ^-^^ - (^Az-^ + 1)1 (A,. + - ^ ■ V WX-z^ + l)(Az-'^ + 7) • 

Let X be such that ^r^r = ^- Then x £ [0, +oo) and by substituting ^rrr for z, we have 

, d\ 2 



a(d;xi, ...,Xrf) < ^ = ad{x). 

+ 1) ^ + 7) ^ (/3A (f±i) % l) ' (a (fd^) % 7 



□ 

Lemma 14. Let (/3,7,A) 6e anti- ferromagnetic, /f (/3,7, A) zs up-to-A unique, then there exists a 
constant a < 1 such that for any integer I < d < A, it holds that aa^x) < a for all x > 0. 

Proof. Fix d to be any positive integer. We characterize the value of x at which ad{x) achieves its 

maximum. We denote that fd{x) = X ^^p^^ . Taking derivative of ad{x) with respect to x, we 
get that 

G'{x) 



a'dix) = d{l - /37) 



(/3/d(x)+l)(/d{x)+7){/3x+l)(x+7) • 



where G{x) = ,sf,(,^+^)(ffJ^l)(B.+^^(r+'y) ^ ^^^ose derivative is 



, fdix) ■ dil - /3j)x f j-/3x^ 7 - /3fdix)' 

G{x)- 



{/3fd{x) + l){fdix) + j){f3x + l)2(x + 7)2 \d{l - M)x {/3fd{x) + l)(/,(x) + 7) 
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As X ranges over [0,oo), the function is strictly decreasing in x and ranges from +00 to 

—00, and the function ?\fl\('f \ \, \ is strictly increasing in x and has a bounded range. Thus, 



the equation 

7-^x2 7 - I3fd{x) 



(4) 



d{l-M)x (/3/rf(x) + l)(/d(x) + 7)' 
has unique solution in (0,oo), denoted by x^- Moreover, it holds that 

> \i Q < X < Xd, 
G'{x)\ = Q ifx = Xrf, (5) 
^< ii X > Xd- 

The same also holds for a'^{x). Thus, for any fixed d, ad{x) achieves its maximum when x = Xd- 
Therefore, for all x > 0, 



adix) < ad{xd) = d{l - p-f) 



Xdfd{Xd) 



{(3xd + l){xd + 7) iPfdixd) + 1) ifdixd) + 7) 

2 



_ o ^ (7 - ^xg) fdjXd) , , . 

^ ''^>- {f3xd + l){xd + j) ' {^-Pfd{xd?)) ^ ^ 

- Oid{Xd)- 

Equation ([6]) is obtained by substituting {Pfd{xd) + 1) {fd{xd) + 7) according to ([1]). 

Let Xd be the positive fixed point of fd{x), that is, Xd = fd{xd)- We then claim that if 7, A) is 
up-to-A unique, then ad{xd) < OLd{xd) for any integral 1 < d < A, To see that this claim is sufficient 

to imply the lemma, note that after substituting Xd = fd{xd), we have otdixd) = ^ (^^^^^^(^^J^-y) = 



And due to Lemma [H if (/3,7, A) is up-to-A unique then there exists a constant c < 1 
such that \fd{xd)\ < c for all integer 1 < c? < A. 

We then prove the claim. Assume that (/3,7,A) is up-to-A unique and 1 < c? < A. It is then 
sufficient to show that ad{x) is decreasing if Xd < Xd and is increasing if Xd> Xd- 
Case 1: Xd < Xd- Due to ®, G'{xd) > 0. Note that 



(/3frf + l)3(xrf + 7)3 \d{l-p-f)xd i/3xd + l){xd + -f) 

Due to the uniqueness, \f'd{xd)\ = (fffl+"i)$r+7) < 1' d(T=kw~, ~ (/3x,+i)(x,+7) > Combining 
with that G'{x) > 0, we have 7 — /3x^ > 0. Since the function fd{x) is monotonically decreasing 
and Xd is its fixed point, 7 — fifdixd)^ > 7 — Pfdixd)^ = 7 — f^x^d — ^- Since Xd satisfies 7 — /3x^ 
and 7 — (3fd{xd)'^ must be simultaneously positive or negative, thus it also holds that 7 — > 0. 

Then both and (^^-i^f^^x)'^) positive and monotonically decreasing in x E [xd,Xd\. 

Therefore, ad{xd) < ad{xd). 

Case 2: Xd > Xd- By the same argument as above, it holds that 7 — [ifdixd)"^ = 7 ~ l^x^d ^ 0' 
7 - f^fdixdf < 0, and 'y - f^xj < 0. Thus both and (^^J^ff^^yi) are negative and 

monotonically decreasing in x G [X(i,Xrf], hence their product is positive and increasing in x G 
[xd,Xd]- Therefore, adixd) < ad{xd)- □ 

Lemma [12] is proved by combining Lemma [13] and [HJ This completes the proof of Theorem [9] 
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Strong spatial mixing on regular trees. 

As a byproduct of our analysis, we prove a strong spatial mixing theorem for regular trees. When 
the graphs G itself is a regular tree. All vertices (except the root) has the same arity. And all ds 
(excerpt the one of the root) that appear in the proof are the same and equal that arity. Then the 
condition that the uniqueness holds on all infinite regular trees of degree up to A can be replaced 
by the uniqueness on infinite A-regular tree. 

Theorem 15. For two-state anti-ferromagnetic spin systems, on any infinite A-regular tree the 
uniqueness implies the strong spatial mixing. 

The same result can be obtained by combining the same theorem for the hardcore model |25j 
and anti-ferromagnetic Ising model [20] and translating the parameters of general two-state anti- 
ferromagnetic spin systems to these models (as discussed in p0ll23]). However, unlike the hardcore 
and the anti-ferromagnetic Ising models, for general two-state anti-ferromagnetic spin systems 
Theorem [15] itself is not sufficient to imply the strong spatial mixing on graphs of maximum degree 
at most A. This is discussed in Section [5] 



4 Algorithmic implications 

In this section we prove Theorem [2] That is, if an anti-ferromagnetic (/3, 7, A) is up-to-A unique 
then there exists an FPTAS for the partition function Z{G) for any graph G of maximum degree 
at most A, and in particular the universal uniqueness implies an FPTAS for arbitrary graph G. 

It is well-known that Z{G) can be computed from p^'^ by the following standard procedure. 
Let vi,. . . ,Vn enumerate the vertices in G. For < i < n, let Uj be the configuration fixing the 
first i vertices vi, . . . ,Vi as follows: (Ti{vj) = (Tj„i(fj) for 1 < j < i — 1 and (Ti{vi) is fixed so that 
Pi = Pr[crj(fj) I o"j_i] > 1/3. In particular, o"„ G {0, 1}^ is a configuration of V . It holds for the 
Gibbs measure of (T„ that /9((T.„) = piP2 ■ ■ ■ Pn as well as that p{o-n) = ^z^g) ' ^^^^ ^i^) — p^i'^"p ' 
where the weight w{an) = I\{u,v)(^E ^rTn{u),an{v)I\v€V b<Tn{v) can be computed precisely for any 
particular it„ in time polynomial in n. Note that pi equals to either Pvl~^ or 1 — Pvl~^ ■ Therefore, 
if p"'^ can be approximated within an additive error e in time polynomial in n and ^, then the 
configurations cjj can be efficiently constructed such that all pi are bounded away from 0, thus the 
product P1P2 ■ ■ - Pn can be approximated within a factor of (1 it ne) in time polynomial in n and 
which implies an FPTAS for Z{G). 



Bounded degree graphs. Let G be a graph whose maximum degree is at most A and v be any 
vertex. A self-avoiding walk tree T = T^p;\^{G,v) can be constructed so that R'q"^ = R'^^. We can 
use the recursive procedure described in Section [3] to compute the upper and lower bounds of R'^^ , 
with the setting that for all the vertices more than t steps away from the root f , the trivial bounds 
< R'^^ < 00 is used. Then the proof of Theorem [9] shows that the recursive procedure returns Rq 
and Ri such that Rq < R'^^ < Ri, and Ri — Rq = 0(a*) for some constant a < 1 assuming that 

(/3,7, A) is up-to-A unique. Note that R^^ = Rq\ = y^^/a • Let po = and pi = ■ Then 

Po < Pv^ < Pi and 
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The recursive procedure runs in time 0(A*) since it only needs to construct the first t levels of 
the self-avoiding walk tree. If A is bounded, by setting t = log„ e, this gives an algorithm which 
approximates pj* within an additive error 0(e) in time polynomial in n and -, which implies an 
FPTAS for Z{G). 

Arbitrary graphs. Let G be an arbitrary graph and v be any vertex. Let T = Tsaw{G,v). We 
use the method of Computationally Efficient Correlation Decay introduced in [TH] to deal with the 
unbounded degrees. Intuitively, using this method we observe correlation decay in a refined metric 
instead of graph distance such that in this new metric a neighborhood of moderate size is sufficient 
to guarantee desirable correlation decay. 

Similarly, we use the recursive procedure described in Section [3] to compute the upper and lower 
bounds of R'^^ , but this time the termination condition relies on a new depth defined as follows. 

Definition 16. Let T be a rooted tree and M > 1 be a constant. For any vertex v in T , define 
the M-based depth of v, denoted Im{v), as such: Im{v) = if v is the root, and Im{v) = ^m{u) + 
\\ogj^{d + 1)] if V is one of the d children of u. 

Let Af > 1 to be fixed. Denote by B[t} the set of all vertices with M-based depth < H. along 
with their children and grandchildren in T. It can be verified by induction that < n^M^. 

The recursion is applied to estimate the B!^ when the current v £ B{i) until v is no longer in B{i) 
in which case the trivial bounds < R"^^ < oo is used. 

Let ey be defined as in Section [3l Repeatedly applying Item [3] of Lemma [TTl without loss of 
generality, we have a path uiU2 ■ ■ ■ in T from the root ui = v to a Uk with lM{uk) ^ ^ and 
^M{uk-i) < (■■, such that e„^. < a[dj] xi, . . . , ) • ^Uj+i for j = 1,2, . . . , k, where dj is the number of 
children of uj and Xi € [0, oo), 1 < i < dj. The key to overcome the explosion caused by unbounded 
degrees is to observe that the contraction a{d;xi, . . . ,Xd) decreases dramatically as the degree d 
grows. 

Lemma 17. Let [(3, 7, A) be anti-ferromagnetic. If (/?, 7, A) satisfies the universal uniqueness, then 
there exist constants a < 1 and M > 1 such that for any integer d > 1, and any Xi G [0, 00), 
l<i<d, it holds that a{d; xi,...,Xd)< ari°SA/(d+i)l . 

Proof. Assume the universal uniqueness of (/3, 7, A). Due to Lemma [T2l there exists a constant a < 1 
such that a{d; xi, . . . , x^) < a. By Item [2] of Lemma [211 the universal uniqueness implies that 7 > 1, 
thus there exists a constant M > 1 such that d ■ ^J^^ < al'iogMC'^+i)! for all d> M. Due to ([3]), it 

holds that a{d; xi,...,Xd) < d- ^ ^ < ari°SA/(<i+i)l fo^ ^h d > M. Note that a{d; xi,. . . ,Xd) < a 

means that a{d;xi, ...,Xd) < a ^'"s a/ ('^+1)1 for d < M. Therefore, a{d;xi, ...,Xd) < a^^°SMid.+i-)] foj. 
all d. □ 

By Lemma [TTl there exist constants a < 1 and M > 1 such that 

k 

i=i 

With the notation used in Section O S is the complement of B{i). Note that all n^'s children are in 
B{£) thus none of them are in S, and by Item [2] of Lemma [21] the universal uniqueness implies that 
7 > 1. Thus by Item [2] of Lemma [TT] it holds that e^j. = 0(1). Therefore, < e^^, • = 0{a^). 
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Let 5i, = Ri — Rq, where Rq and Ri are the bounds returned by the recursive procedure such 
that Ro < R^'^ < Ri. By the same analysis as in Section [3l 5^ = tt^: = 0{ey) = O(a^). Then 

by ([7]), the marginal probability is approximated within an additive error of 0{a^). The running 
time of the recursion is 0{nB{i)) = 0{n^M^). By setting £ = log^^e, we have an algorithm which 
approximates within an additive error of 0(e) in time polynomial in n and ^, which implies an 
FPTAS for Z{G) for arbitrary graph G. 

Heterogeneous spin systems. Our analysis in last and this sections actually holds for hetero- 
geneous spin systems which allow that each vertex v has a distinct constant external field > 0. 

Theorem 18. For a two-state anti-ferromagnetic heterogeneous spin system with parameters j3, 7, 
and external field at each vertex v, for any finite A > 2 or A = 00, if for all v the (/3,7, A^,) 
is up-to-A unique then the spin system is of strong spatial mixing and has FPTAS on graphs of 
maximum degree at most A. 

5 Non-monotonicity of general two-state spin system 

In this section, we prove the following theorem. 

Theorem 19. There exist two-state anti-ferromagnetic spin systems which exhibit strong spatial 
mixing on infinite d-regular tree but does not exhibit strong spatial mixing on all graphs of maximum 
degree at most d. 

In a seminal work [25], Weitz proved that for the hardcore model the strong spatial mixing on 
an infinite d-regular tree implies the strong spatial mixing on graphs of maximum degree at most d 
(Theorem 2.3 in |25j). He further remarked that the same implication holds for all two-state spin 
systems. That is, 

for any two-state spin system, strong spatial mixing on an infinite d-regular tree implies 
the strong spatial mixing on graphs of maximum degree at most d. 

This claim played important roles in current understanding of correlation decay in two-state spin 
systems as well as devising FPTAS for such systems. An algorithmic form of this claim was cited 
in [20] as a theorem for all two-state spin systems (Theorem 2.8 in [20]) and was proved for the 
anti-ferromagnetic Ising model. It was raised as a conjecture in [21] whether this claim holds for 
multi-state spin systems. 

Here we clarify that this claim holds only for the two-state spin systems under limited settings 
but does not hold for all general two-state spin systems. This disproves the conjecture in [21j and 
shows that the common belief that the d-regular tree represents the worst case for strong spatial 
mixing among all graphs of maximum degree at most d cannot be generalized to general two-state 
spin systems. We first describe a region that the claim is true. 

Lemma 20. For < /3,7 < 1, the strong spatial mixing on infinite d-regular tree implies the strong 
spatial mixing on trees of maximum degree at most d. 

Proof. Given a rooted tree of maximum degree at most d, for each vertex of k children with k < d— 1, 
we can attach d — 1 — A: dummy children with fixed (distributions of) spin states. This is the method 
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used in [25] and [20]: for the hardcore model the dummy children are fixed to be unoccupied and 
for the anti-ferromagnetic Ising model the dummy children are of uniform distributions over spin 
states. In both cases, the dummy children have no effect on their parent. In general, we fix the 
distribution to be ipo,pi) at each dummy child satisfying 

PQ+pi = 1, 
f3po+Pi =Po + IPi- 

When < /3,7 < 1, this system has solutions in < pojPi ^ 1- With the ratio Ri at the i-th child, 
1 < i < /c, and Ri = po/pi for the dummy children k < i < d — 1, the ratio at the parent is given 
by the recursion 

which is identical to the original quantity. □ 

Due to the self-avoiding walk tree construction (Theorem [6]), it holds that for < /3,7 < 1, 
strong spatial mixing on infinite d-regular tree implies the strong spatial mixing and the FPTAS 
on graphs of maximum degree at most d. 

For < /3,7 < 1, the spin system shows the following monotone property: the uniqueness 
on infinite d-regular tree implies the uniqueness on all infinite regular trees of smaller degree. 
This can be verified by the following reasoning: due to Theorem 1151 on infinite d-regular tree the 
uniqueness implies the strong spatial mixing, which for < /3, 7 < 1, implies the strong spatial 
mixing (including the uniqueness) on all infinite regular trees of smaller degree. 

There exist two-state anti-ferromagnetic spin systems which are non-monotone. We can choose 
anti-ferromagnetic (/3,7,A) satisfying that 7 > 1 and A > Ac(/3,7), where Ac(/5, 7) is the critical 
threshold for universal uniqueness given in Item [8] of Lemma [2TJ According to Item [3] of Lemma 
[211 since 7 > 1 the uniqueness holds on d-regular trees for all sufficiently large d. On the other 
hand, according to Item [8] of Lemma [211 since A > Ac(/S,7), there exists a finite d' such that the 
system is non-unique on d'-regular tree. 

For such non-monotone systems, due to Theorem 1151 the uniqueness implies the strong spatial 
mixing on d-regular tree for sufficiently large d, but the strong spatial mixing does not hold on a 
regular tree of smaller degree (because of the non- uniqueness on the smaller tree). Therefore the 
implication between the strong spatial mixing on d-regular tree and on graphs of maximum degree 
at most d does not hold for general two-state spin systems. 

Acknowledgment We would like to thank Alistair Sinclair and Piyush Srivastava for the helpful 
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A The Uniqueness Thresholds 

The following lemma translates the uniqueness condition into its various threshold forms. The 
lemma (except for item [2]) is not used in the proofs of the main results but is used in the inter- 
pretation of the main results and comparisons to the previous results which are mostly stated in 
threshold forms. 

Lemma 21. Let (/3,7,A) be anti- ferromagnetic. 

1. (/3,7,A) is up-to-2 unique. 

2. If^^^, then the uniqueness does not hold on infinite d-regular tree for all sufficiently large d. 

3. If J > 1, then the uniqueness holds on infinite d-regular tree for all sufficiently large d. 



18 



4- For any A, there exists a critical threshold -fc = idf^, ^) such that (/3,7,A) is up-to-A 
unique if and only if j G {"fc, In particular, 7c(/3, A, oo) > 1 and jdf^, A, oo) = 7c(/3, A, A) 
for some finite A. 

5. If 13 = 0, for any A, there exists a critical threshold Ac = Ac(7, A) = mini<(i<A (J_i)d+i such 
that (/3,7, A) is up-to-A unique if and only if X £ (0, Ac). 

6. If ^/JJJ > then for any external field X, (/3,7,A) is up-to-A unique. 

7. If (3 > 0, for any A that a/37 < there exist two critical thresholds Ac = Ac(/9,7, A) and 
Xc = Ac(/3,7,A) such that (/3,7,A) is up-to-A unique if and only if X £ (0, Ac) U (Ac,oo). In 
particular, when 13 = ^, it holds that Xc - Xc = 1; thus {f3, f3, A) is up-to-A unique if and only 
^/|logA| > log Ac. 

8. If 13 > and 7 > 1, there exists a absolute positive constant Xc = Xc{f3,^) such that (/3,7,A) 
is universally unique if and only if X £ (0, Ac). 

^xa+-y ) ~ fdi^d) be the positive fixed point of /d(x). Then 

l^'^^'^'^^l- (/3x, + l)(x, + 7)- 
Let (/3,7, A) be anti-ferromagnetic. That is, < /3 < 7, 7 > 0, and /37 < 1, thus /3 < 1. 

1. It is easy to verify that for /?7 < 1, (/5a; + l)(x + 7) — (1 — f3^)x > for any positive x > 0. 

(1-/37)^1 

(/3xi+l)(xi+7) 



Therefore, when d = 1, we have |/i(xi)| = /o-^V fj-^Y x < 1, which means that (/3,7,A) is 



up-to-2 unique. 
2. For ah sufficiently large d, it holds that 



^pd I d \ (i(l-^7)-3 , 

'^^^P ^ (1-/37)^-3 ) < / ' 

Aexp —-77- > 



d{l-pj)-3j d{l-(3j)-^' 
By contradiction, suppose that |/d(xd)| = < 1- Then, 

d{l-(3j)xci _ d{l-(3j) ^ d(l-/37) 



{pxd + l){xd + 7) pxd + g + {l + ^j) - pXd + g + 2- 



Case.l: x^ > 7. Then ^ <1. Thus, 



- Pxd + ^^+2- (3x^ + 3 

-) 

/3 



which implies that x^ > ^(L_^l)_^, However, it holds that 



Xd + 1 J \ Xd J V d{l- /3'y) 



(l-|97)<i-37 /3 



a contradiction. 
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Case. 2: Xd < 7. Then (5xd < (3j < 1. Thus, 

d(l-/37) ^ rf(l-/37) 

which imphes that < . However, it holds that 



Xd + l ) - {xd + lY - \ 41-/57) -3 

^ \ ( ^ ^ ^ 

> A • exp — — — > 



d{l-pj)-3j 41-/37) -3' 
a contradiction. 

3. Let 7 > 1. The fixed point = A ^ ^ ^ S'l^d 

which is strictly less than 1 for all sufficiently large d. Thus the uniqueness holds on infinite 
(i-regular tree for all sufficiently large d. 

4. We first show that there exists a critical threshold 7c = 7c(/3, A, A) such that (/?, 7, A) is up- 
to-A unique if and only if 7 G (7c, ^). It is sufficient to show that if an anti- ferromagnetic 
{(3, 7, A) is up-to-A unique then 7', A) is up-to-A unique for any 7' > 7 and 13^ < 1. 

Recall that Xd is the positive fixed point of fd{x) = A ^^^p^^ . Also let x'^ denote the positive 

solution to X = A ' 



We first show that x'^ < Xd- By contradiction, assume that x'^ > Xd- Since for anti- 
ferromagnetic (/3,7,A), fdix) is monotonically decreasing, we have 



Xd = X{ ^ " > A > A 

^ Xd + 7 J \Xd + l J \x'd + J 



a contradiction. 
Since 



V x'^ + i J \ xd + i 



we have 

(i-/3y) ^ (1-/37) 

For < Xrf, it also holds that 



/3x;^ + i /3+1 p + j- pxd + i 
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Multiplying the above two inequalities together, we have 



i(3x'^ + l){xd + 7') + 1) (xrf + 7) ■ 

Note that these are the absolute derivatives at the respective fixed points when the parameters 
are (/3,7, A) and (/3,7',A). Therefore if (/3,7,A) is up-to-A unique then (/3,7',A) is up-to-A 
unique. 

Due to Part[5]of the lemma, if 7 < 1, \f'^{xd)\ > 1 for all sufficiently large d, thus 7c(/3, A, 00) > 
1. And due to Part [3] of the lemma, for any 7 > 7c(/5, A,oo) > 1, \ fdi^d)\ is arbitrarily close 
to as d grows to infinity, thus 7c(/3, A, 00) = 7c(/5, A, A) for a finite A. 

5. When /3 = 0, = (^g^^^^gf^^) = £f^, the uniqueness condition < 1 is 

equivalent to that x^ < (here we assume d > \ since due to Part [T] of the lemma, for 

d = 1 the uniqueness always holds). Recall that = A i^-^r^j^ ■ Then < 1 if and 

only if 

A = Xd{xd + 7) < z^^wTT- 



(d-l) 

Let Ac = Ac(7, A) = mini<^<A (J_i)f+i • It holds that (/3, 7, A) is up-to-A unique if and only 
if A G (0,A 



6. We note that \fd{xd)\ = (/j^^+i'^gj+.y) is not monotone in Xd- It achieves its maximum value 
at Xd = yj^. Therefore, if for any 1 < d < A, (^l+7f(I+7) < 1 for x = yj^, then (/3,7, A) is 
up-to-A unique for any A. This condition holds when > for all 1 < d < A, i.e. when 

7. The following lemma is needed in the proofs of this part and next part. 
Lemma 22. Given that < /3 < 7, /37 < 1, and \fWl ^ ^TT' define 



-1 - /37 + d{l - /37) - Vi-l - /37 + ^1 - - 
Md) = , 



-1 - /37 + d{l - + - /37 + ^(1 - h)? - 4/37 

X2{d) = . 

which are the two positive roots of equation (^'f'^7j(x+7) ~ 



Let \i{d) = Xi{d) (|,r(3)+T 



Xt(d)+7 



d 



1,2. 



(a) xi{d) is monotonically decreasing in d; X2{d) is monotonically increasing in d. 
(h) \2{d) is monotonically increasing in d and goes to infinity as d grows, 
(c) If J > 1, Xi{d) achieves its minimum value at a finite d. 
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Proof. Note that {-1 - l3j + d{l - - 4/37 > if V^f < which means that xi{d) 

and X2{d) are weh-defined. 

(a) It is obvious that X2{d) is increasing in d and is unbounded as d grows. Note xi {d)x2{d) = 
^, so xi{d) is decreasing in d and tends to as (i grows to infinity. 

(b) For < /? < 7, /37 < 1 and y/JJ^ < we have /3 < y/^ < ^ < 1 and is 
increasing in x. Moreover, 

X2(c?)+7 ^ 1 (rf - 1)(1 - /37) + v/(-l - /37 + ^(1 - /37))' - 4/^7 

/32;2((f) + l /? ' (d + 1)(1 - /37) + v^(-l - /37 + d{l - p-f)f - 4/37 

^ d+1 (d-l)(l-/37) _^ 

- d-1 (d + l)(l-/37) 

Since X2(d) is also increasing in d and unbounded as d grows, we can conclude that 
= X2id) {fl^y is increasing in d and is unbounded as d grows. 

(c) It is sufficient to prove that Xi{d) is unbounded as d grows. When 7 > 1, we have 



Xi{d) = xi{d) 



xi{d) + 7 



(3xi{d) + l 
> x,{d)-Y 



-1 - /37 + d{l - /37) + V(-l-/37 + ^i(l-/37))'-4/37 
which is obviously unbounded for large d. 

□ 

Recall that < Let A be the smallest integer d satisfying that -v/^ < Then 

we have -y/^ > 3^ for all 1 < d < A and < for all A < d < A. According to Part 
[6]of this lemma, for all A > 0, (/3,7,A) is up-to-A unique, i.e. < 1 for 1 < d < A. 

It remains to analyze the for such d that A < (i < A, i.e. \/^ < By Lemma 

d{l-(3-y)x 
(/3x+l)(a;+7) 



the equation w"^|f ^ = 1 has two positive roots xi{d) and X2(d) as given in Lemma 



It then holds for such d that = < 1 if and only if x^ < xi{d) or 

Xd > X2{d). Note that x i^ ^^^^ is monotonically increasing in x for any fixed d, thus 

Therefore, < 1 if and only if A < Ai(d) or A > \2{d). 
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For any A that < let 

Ac = Ac(/3,7, A) = _min Xi{d), 

A<d<A 

Ac = Ac(/3,7,A) =_max A2(d). 

A<ci<A 

It holds that (/3, 7, A) is up-to-A unique if and only if A G (0, Ac) U (Ac, 00). 
We then show that when /3 = 7, it holds that Ac-Ac = 1. First, it is easy to see that when /3 = 7, 
xi{d)x2id) = ^ = 1, thus Xi{d) ■ X2{d) = xi{d)x2{d) {j^^^ ■ plf{d)+i ) = 1 because /3 = 7 
and xi{d)x2{d) = 1. Due to Part [7bl of Lemma [22l \2{d) is monotonically increasing thus 
Xi{d) = l/\2{d) is monotonically decreasing. Therefore, Ac = ™-^'^~^<:d<cA 

X2{d) = A2(A-1) 

and Ac = ™iii5-<^<^ Ai(d) = Ai(A — 1) and it holds that Ac • Ac = Ai(A — 1) • A2(A — 1) = 1. 
Then 7, A) is up-to-A unique if and only if A G (0, 1/Ac) U (Ac, 00), i.e. | log A| > log Ac. 

8. Let A be the smallest integer d satisfying that \fWy ^ ^^pr- Then by Part [71 for any A, 
(/3,7, A) is up-to-A unique if and only if A G (0, Ac) U (Ac, 00), where Ac = min^<^^^ Xi{d) 
and Ac = max-^^^^^ A2((i). By Part [7a] and Part I7bl of Lemma [2^ for fixed /3,7, \2{d) goes 
to infinity as d grows. Therefore, Ac is unbounded for A = 00. 

If 7 > 1, by Part[7c|of Lemma [22l there is an absolute positive constant Ac = min^^,^ Xi{d) 
such that (/3, 7, A) is universally unique if and only if A G (0, Ac) U (Ac, 00) = (0, Ac). 

□ 

B Proof of Lemma [11] (a Mean Value Theorem approach) 

We prove Lemma [TTJ The notations are defined in Section [3l 

1. Due to the Mean Value Theorem, there exists an i? G [Ri,,Rv + 8^] such that 

= ^{R^ + 6^) - ip{Ri,) = ip'{R) • 6^, = ^{R) ■ 6^. 

2. Suppose that each vertex has at most k children. Then due to the assumption either k is 
bounded or 7 > 1. 

If V A, then (5^, < Rv + = f{Ry, , . . . , i?„J < /(O, . . . , 0) = ^, where < d < A;. If 7 > 1, 



then R^ + 6v < X = 0(1); and if k is finite, then R^ + 6^ < max{;^, A} = 0(1). 
Due to the Mean Value Theorem, there exist Ri G [Rvi,Rvi + (^nj, 1 < i < d, such that 
ev = ip{f{Rvi,. . . ,-Rt.J) -ip{f{Rv^ +5v^,...,Rv^ + 6^^)) 

= -V(^ (^f{Ri,...,Ra)y{6,„...,6,,) 

1 
2 



(1-/37)- (/(i?i,...,i?.)j^ _ ^ 5., 



d 



< 



/ /(0,...,0) y.^,. 
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If Vi S for all Vi, then due to the above argument, S^. < ^ for free Vi, where < di < k is the 

number of children of vf, and trivially 6^. = for fixed Vi. Therefore, < \J Ylt=i yJ^fi — 

\2d'y~^~ -max ^-^^ l|- If 7 > Ij then it is easy to see that d'y~~^ is bounded by a constant, 
thus it holds that = 0(1); if k is bounded, then d < A; is also bounded, thus clearly 
e. = 0(1). 

3. We then analyze the stepwise contraction of e„. Define that y„ = (p{Rv) and accordingly 
y^. = ifiRv,), 1 < i < d. Then y„ + e„ = ip{Ry + 5„) and y^,. + = (/^(i?^^ + 1 < i < d. 
We have 

Vv + e^ = 93(/(93"^(y„J, ...,93"^(y„J)). 

Apply the Mean Value Theorem. There exist yi € [yviiUvi + CdJ and corresponding Ri G 
[i?,;., + 6vi] satisfying yi = ip{Ri), I < i < d, such that 



-c^(/(i2;,...,ir,))-^^-L-.6, 



2 



(f,xuU^ + 1)' (Ant, fJ7 +7)' <^^- + 

(i-w(Ant.fg)' f Hi 

< max {e„^ | j ' ■ } 



^--"^ (^/3Anti^^ + l)' i=i (/^^^ + l)'(^^ + ^)• 



a(d;i^l, . . . • max{e„.}. 

l<i<(i 



C Finding a good potential function heuristically 

Perhaps the most mysterious step in our proof is the choice of the potential function <I>(x) = 

. ^ As in many cases where potential analysis is applied, there is no standard routine 

-^x(/3x+l) (x+'y) 

for searching for a suitable potential function. On the other hand, it is quite unlikely that we 
can just guess such a fairly complicated formula without any hints. Here we present a heuristic 
approach which leads us to the discovery of a good potential function. This part is not rigorous 
and logically unnecessary for the soundness of our result. Nevertheless, this heuristic approach is 
general and interesting enough and may find its applications in other scenarios, thus deserves an 
exposition here. 

The heuristics consists of three steps: 

1. Find a necessary condition for the potential function, which is an equation related to the 
potential function at one point. 
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2. (heuristic step) Enhance the condition by assuming that the equation holds for the whole 
range of the variable, which gives us a differential equation. 

3. Solve the differential equation and get a potential function. 

We first assume that the system is at the boundary of uniqueness. This means that d is the 

and X = f{x) is the positive fixed point. 



critical degree and f'{x) = —1, where f{x) = A f^^pp 



We have the following identities: 



X = X 



/3x + 1 



and 1 = Ad(l - /Sj) 



{f3x + l 



X + 7 y ' ' " (x + 7) 

which together implies another identity 

d(l-/37)x = + l)(x + 7). 



d+l 



(8) 



The goal is to find a potential function such that a{x) = \ f'ix)\^^^y- < 1 for all x. On the 

1. So a{x) achieves its maximum when 



1 • 



<i.(i) 



other hand, we have a{x) = \f'{x)\-^^ - j- • -^py 
X = x. As a differentiable function, it holds that a'{x) = 0, that is 



fix) 



$(x) 



0. 



We have the following calculation: 



^ [/'(x)$(/(x))]'$(x) = /'(x)$(/(x))$'(x)|^^, 

^ [f"{xMf{x)) + /(x)$'(/(x))/'(x)] $(x) = /'(x)$(/(x))$'(x) 

/"(x) <I>'(x) 



where we use the fact that x = f{x) and /'(x) = —1. 
Taking the second derivative of /(x), we have 

(ffr 4- 1)^-'^ 

fix) = \d{M - 1) ^ty)i+2 - + 7) - + + 1)) • 

This equation already gives an identity for But it is too complicated. We may further simplify 
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the above formula at the point x = x. 



f"{x) = \d{f^^-i)^-^±^[{d-mx+^)-{d+mx+i)) 

^ ((d + l)(/3x + l)-((i-l)/3(x + 7)) 



(/3x + l)(i + 7) 
d+l id-l)P 
£ + 7 /3x + 1 

d{l - ^ _}_ ^ /3 



(/3x + !)(£ + 7) £ + 7 px + l 

1 1 /3 
+ + 



So we have 



X X + 7 /3x + 1 

We apply the heuristics and assume that Q simply holds for all x. This gives us a differential 
equation 

(ln($(a:))y = -i(i + -l-+ ^ ^ 



2x x + 7 fix + 1 
The solution of this differential equation is 

ln($(x)) = ln(x(3; + 7)(/3x + 1)) + Ci, 

which gives that 



V^x(;Sx + l)(x + 7)' 

where Ci, C2 are some constants. This gives the potential function we used in the paper. 

There are also other equations which hold for the fixed point x. Choosing which equation for 
X to heuristically extend to all x may affect the potential function we obtained. For example, ([8]) 
implies that 

1 _ 1-/37 _ 1 P 
dx {(3x + l){x + j) x + 7 /3x + l' 

thus we can rewrite /"(x) as 



This gives that 



X x + 7 /3x + 1 dx I3x + 1 
(ln($(x)))'= ^ 



2dx /3x + 1 

We can also treat this equation as a differential equation for variable x and solving it gives us 

^{x) = -, 

X 2d (/3x + 1) 

which is the potential function used in |15| . 
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